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1 . 0 Introduction 


There  are  two  claims  about  meteorological  satelli  tes  which  are  not 
likely  to  be  challenged.  The  first  claim  is  that  they  have  provided  us 
with  an  abundance  of  data,  and  the  second  is  that  the  data  have  not  been 
converted  into  information  in  a very  efficient  manner.  The  reasons  for 
this  inefficiency  are: 

1.  Hard-copy  imagery,  although  a valuable  tool  in  the  hands  of  a 
skilled  analyst,  suffers  from  inconsistency  in  interpretation,  is  expensive 
to  distribute  widely,  and  is  time-consuming  to  process  in  the  variety  of 
formats  the  analyst  might  need. 

2.  Techniques  and  the  hardware  to  process  data  objectively  are  in- 
adequate to  handle  the  overwhelming  volume  of  data  transmitted  by  both 
polar  orbiting  and  geosynchronous  satellites. 

The  difficulties  with  hard-copy  imagery  are,  to  a large  extent,  being 
overcome  with  interactive  systems  such  as  McIDAS  (Man-ccnputer  Interactive 
Data  Access  System)  which  combine  the  speed  of  electronic  image  processing 
with  the  judgement  of  the  skilled  analyst. 

Objective  processing,  on  a rudimentary  level,  has  been  taking  place 
at  the  Air  Force  Global  Weather  Central  (AFGWC)  since  1971  (Cobum  1971)  to 
support  the  specification  and  forecasting  of  global  cloud  cover.  Satellite 
data  are  combined  with  surface,  radiosonde  and  aircraft  observations  to  pro- 
duce a 15-level  global  cloud  distribution  with  a horizontal  resolution  of 
25  miles  square.  This  product,  named  3-D  Nephanalysis,  is  used,  ancng 
other  ways,  as  the  initial  condition  for  the  cloud  forecast  in  the  global 
miner ical  forecast.  The  AFGflC  is  the  only  operational  group  that  atbanpts 
to  exploit  satellite  imagery  in  an  automated  process. 


Despite  AFCJ'C's  leadership  in  the  application  of  satellite  data,  re- 
quirements placed  on  them  by  their  customers  exceed  their  processing  capa- 
bility. Presently,  anoothed  mode  data  at  approximately  3 nm  resolution 
are  used  as  the  satellite  input  to  the  3-D  Nephanalysis  although  fine  mode 
data  at  1/3  rm  resolution  are  sensed  by  the  satellite.  The  fine  mode  data 
are  displayed  as  hard  copy  but  are  not  otherwise  used  because  until  recently 
the  AFOC  computers  were  completely  inadequate  to  operate  in  real  time  with 
the  rates  and  volume  of  data  involved. 

This  report  discusses  seme  initial  efforts  the  Regis  College  Research 
Center  directed  toward  the  development  of  methods  of  improving  the  3-D 
Nephanalysis  by  the  use  of  fine  mode  data  and  eventually  the  development 
of  a model  that  uses  only  satellite  data  to  define  the  cloud  cover. 
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2.0  Approach 


2.1  The  Nature  of  the  Problem 

Inprovanent  of  operational  cloud  specification  at  AFQWC  by  means  of 
the  utilization  of  fine  mode  data  poses  two  problems.  The  first  problan  is 
hew  to  handle  the  enormous  quantities  of  data  involved.  The  fine  mode 
sensor  of  the  Block  5-D  system  generates  over  100  times  as  much  data  as  the 
high  resolution  sensor  in  the  Block  5-C  system. 

Figure  1 illustrates  what  this  100-fold  increase  in  data  means  in 
terms  of  the  amount  and  quality  of  the  data.  The  images,  photographed  on 
the  Mo  IDAS  screen,  are  Block  5-C,  fine  mode  data  processed  so  as  to  simulate 
Block  5-D  anoothed  mode  (3  nm)  data  in  Fig.  la  ard  fine  mode  (1/3  nm)  data 
in  Fig.  lb.  Each  of  these  displays  contain  149,500  6-bit  picture  elements. 
The  ratio  of  the  area  covered  by  these  2 million  bits  of  fine  mode  data  is 
illustrated  by  the  size  of  the  outlined  box  to  the  entire  display  in  Fig.  la. 

The  second  problem  of  providing  improved  information  content  using  sig- 
nificantly larger  amounts  of  data  can  be  approached  by  developing: 

a.  better  software 

b.  better/more  hardware 

c.  a combination  of  improved  software  and  hardware 

Considering  the  magnitude  of  the  data  involved,  sophisticated  programming 
alone  could  not  cone  close  to  alleviating  the  problem.  More  and  larger 
computers  are  an  expensive  way  to  go.  Ihe  approach  selected  by  ATOWC  and 
AFGL  was  to  utilize  sophisticated  analytical  techniques  which  are  adaptable 
to  hard  wired  computational  devices.  Specifically,  this  report  discusses 
preliminary  attempts  to  utilize  coefficients  of  Fast  Fourier  Transforms  (FFT) 
as  indices  of  cloud  types.  Although  any  success  in  this  approach  may  not  be 
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Fig.  la  Block  5-C  fine  mode  (1/3  ran)  data  processed  so  as  to  simulate 
Block  5-D  smoothed  mode  (3  ran)  data.  The  outlined  area  is 
over  Centred.  America.  The  photograph  was  taken  from  the  AGFL 
Me IDAS  display. 


Fig.  lb  Block  5-C  fine  node  (1/3  nm)  data  of  the  area  outlined  in  Fig. 

la.  Note  the  increase  in  detail  between  this  photograph  and 
the  one  in  Fig.  la.  The  photograph  was  taken  from  the  APGL 
McIDAS  display. 
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ironed iately  applicable  to  AFGWC  operations,  it  could  influence  a decision 
to  invest  in  "black  bonces"  to  perform  the  enormous  amount  of  "number  crunch- 
ing" involved  in  real-time  FFT  computations . 

2.2  The  3-D  Nephanalysis  Concept 

^he  3-D  Nephanalysis  model  new  being  used,  cloud  characteristics 
are  defined  in  boxes  approximately  25  miles  square  over  the  earth.  All 
available  data  - satellite,  surface,  radiosonde  and  aircraft  - are  used  as 
input.  Hie  output,  insofar  as  the  input  allows,  specifies  the  amount  of 
cloudiness  at  15  levels  in  the  atmosphere ; the  type  of  cloud  identified  as 
cunulus,  stratus,  cumulonimbus,  etc;  cloud  base  and  top;  and  type  of  precip- 
itation. Over  remote  areas  where  only  satellite  data  are  available,  only 
the  amount,  type,  and  tops  of  the  clouds  can  be  specified.  When  contradic- 
ting observations,  such  as  a difference  in  the  cloud  amount  reported  by  a 
ground  observer  and  the  satellite  algorithm,  are  received,  the  model  makes 
a decision  based  on  a data-reliability  decision  tree. 

2.3  Scope  of  Investigation 

This  investigation  is  concerned  with  just  that  part  of  the  3-D 
Nephanalysis  that  treats  satellite  data.  It  is  called  the  satellite  pro- 
cessor and  is  the  software  that  converts  the  satellite-observed  brightness 
and/or  temperature  into  cloud  categories.  In  this  report  we  will  be  dealing 
only  with  Block  5-C  visual  data  because  of  seme  unanticipated  delays  in  the 
availability  of  Block  5-D  data. 


3.0  The  Satellite  Processor 
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The  satellite  processor  stage  of  the  3-D  Nephanalysis  currently  doc- 
umented by  AFS4C*'1 2  uses  both  visual  and  infrared  data.  The  visual  data, 
which  we  are  using  exclusively  at  this  time,  generate  two  cloud-identifica- 
tion  indices.  The  first  index,  V,  is  the  average  deviation  of  the  64 
brightness  values  in  each  25  mi.  sq.  box.  The  second  index,  A,  is  an 
adjusted  average  brightness: 

A = (B  - M7)Fb  , 

where: 

B is  the  average  brightness  of  a 25  x 25  rm  box, 

M7  is  a seven-day  average  of  background  brightness, 

Fg  is  an  empirically  derived  background  brightness  range  factor. 

Figure  2a  is  a graph  of  Fg  as  a function  of  background  brightness.  Figure  2b 
is  a graph  of  the  adjusted  average  brightness  (A)  incorporating  the  back- 
ground brightness  range  factors  (Ffi)  versus  cloud  cover  for  the  ESSA-9 
satellite.* 

In  the  visual  data  processor  the  type  of  cloud  is  determined  from  the 
adjusted  average  brightness  (A)  and  the  average  deviation  (V) , using  the 
following  set  of  conditions: 

1)  If  V = 3,  broken  cloud  conditions  or  variable  cloud  tops  with  the 
type  defined  as  Cumulus. 

2)  If  V < 3 and  A > 56,  the  type  of  clcud  is  defined  as  Cumulonimbus. 

1.  Maj.  Coburn,  A.R. , "Improved  Three  Dimensional  Nephanalysis  Model", 

APCTC  Technical  Memorandum  71-2,  1 June  1971. 

2.  Capt.  Fye,  F.K. , and  Capt.  G.L.  Logan,  "AFGWC  Satellite  Based  Analysis 
and  Prediction  Programs",  Proc.  7th  Technical  Eacl.  Conference, 
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3)  If  V > 3 and  A < 56,  the  type  of  cloud  is  defined  as  Stratus. 
Infrared  data  are  not  classified  into  cloud  types  but  are  used  to  provide 
cloud-top  height  input  for  the  3-D  Nephanalysis.  Information  that  is  given 
by  both  the  visual  and  infrared  data  processors  are  put  into  a data  has**  for 
subsequent  use  by  3-D  Nephanalysis . 

It  is  possible  to  obtain  the  same  infomation  using  FFT  output  pro- 
cessing. With  the  present  data  flow  structures  at  AFGWC,  incorporation  of 
FFT  output  processing  will  not  require  major  data  flow  alterations.  From  a 
software  standpoint,  incorporation  of  FFT  output  processing  can  be  achieved 
by  adding  a subroutine  to  the  present  satellite  processor.  Depending  upon 
the  manner  in  which  the  FFT  output  is  processed,  it  is  possible  to  extract 
from  the  imagery  even  more  information  than  the  current  visual  and  infrared 
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processors  give.  As  demonstrated  by  Booth  , FFT  output  processing  of  dual 
channel  (IR  and  visual)  imagery  can  yield  cloud  classification  for  a larger 
number  of  cloud  types  than  is  new  available  at  APGMC.  Operating  on  the  FFT 
output  in  a different  fashion  will  extract  information  about  the  alignment 
of  cloud  bands.  If,  however,  the  additional  information  is  to  be  incorpor- 
ated, many  data  flew  structures  within  3-D  Nephanalysis  ani  the  satellite 
processor  will  require  modification. 


The  Fourier  Transform  in  the  past  has  been  utilized  for  characterizing 
linear  systems  and  for  identifying  the  frequency  components  making  up  a 
continuous  vraveform.  However,  when  the  waveform  is  sampled,  or  the  system 
is  to  be  analyzed  on  a digital  computer,  it  is  the  finite,  discrete  version 
of  the  Fourier  Transform  that  must  be  used.  The  Fast  Fourier  Transform 
(FFT) , as  developed  by  Cooley  and  Tukey,  is  simply  an  efficient  method  for 
computing  the  discrete  Fourier  Transform. 

The  Fourier  Transform  for  a continuous  waveform  can  be  written  as: 

la)  Y (f ) = X(t)exp(-2irift)dt 

lb)  X(t)  = /l  Y(f)exp(+2mift)df 

where 

Y(f)  is  the  frequency  domain  function, 

X(t)  is  the  time  domain  function, 
i = / -1 

The  analogous  discrete  Fourier  Transform  pair  can  be  written  as: 

1 N-l 

2a)  Y(j)  = - j^q  X (k) exp ( -2tt  jkiytJ) 

2b)  X(k)  = V Y(j)exp(+2irjki/K) 
j=0 

where,  Y(j)  is  analogous  to  Y(f)  in  la,  X(k)  is  analogous  to  X(t)  in  lb, 
and  N is  the  nurber  of  points  sanpled.  The  discrete  Fourier  Transform  can 
be  computed  from  equations  2a  and  2b  directly,  but  it  is  slower  than  using 
the  FFT.  Par  an  N point  transformation,  where  N=2m  (m  is  a positive 
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integer) , the  number  of  complex  additions  and  subtractions  is  given  by  N . 


For  the  FFT  the  lumber  of  complex  additions  and  subtractions  is  reduced  to 


(N/2)  log^-  Pear  N = 1024  this  represents  a servings  of  over  200  to  1. 

The  output  of  an  FFT  is  simply  the  coefficients  of  a Fourier  series 
representation  of  a function.  A one-dimensional  Fourier  series  representa- 
tion of  a function  reduces  to  a simple  harmonic  expressed  by: 


2niTx. 


3)  F(x)  = Aq  + E^  (AjCost^p)  + BjSi 


Sin( 


2nmc 


■)  ) 


where. 


T = period  of  the  function 


x = sampled  domain 


j = phase 

and  AQ,  A^ , and  are  the  coefficients.  For  a one-dimensional  FFT  of  an 

8-pixel  scan  line  the  output  consists  of  the  .17  coefficients  Aq,  Aj,  and  Bj , 

with  j ranging  from  1 to  8.  In  a two-dimensional  FFT  of  an  8 x 8 pixel  area 

the  output  would  consist  of  the  129  coefficients  A^  h.  . , and  B.  . with 

j*1  3,i 

j and  i ranging  from  1 to  8.  The  two-dimensional  FFT  coefficients  are  then 
used  to  oenpute  the  power  spectrin.  The  proper  processing  of  spectra  will 
extract  information  from  the  imagery  field. 


5.0  Normalized  Average  Power  Spectra  Computation 

When  a power  spectrum  is  confuted  from  2-dimensional  Fourier  Transform 

coefficients,  a new  set  of  parameters  is  generated  in  response  to  the  spatial 

distribution  and  frequency  characteristics  of  the  data  field.  Both  Sikula4 
2 

and  Booth  demonstrated  that  some  cloud  types  generate  reproducible  power 
spectra  parameters,  i.e. , identifiable  spectral  signatures.  These  spectral 
signatures  indicate  size,  shape  and  texture  of  a particular  cloud  type.  The 
power  spectrum  can  be  interpreted  as  the  amount  of  variance  in  an  imagery 
field  due  to  clouds  of  various  sizes.  Booth  demonstrated  that  by  proper 
utilization  of  power  spectra  one  can  classify  clouds  by  using  a discriminant 
function  once  the  boundaries  of  the  cloud  classes  were  determined. 


As  an  initial  step,  power  spectra  were  confuted  by  annularly  integrating 
the  FFT  coefficients.  It  is  advantageous  to  ccrpute  power  spectra  in  annular 
rings  since  the  FFT  retains  the  property  of  a periodic  function  allowing 
duplication  of  FFT  output  as  a half-range  expansion.  Figure  3 schematically 
illustrates  this  concept.  Each  two  digit  group  seen  in  Figure  3 represents 
a two-dimensional  FFT  coefficient  of  power.  Fewer  is  defined  as 
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if  A.  . is  the  real  part  (Cosine  term  coefficient)  and  B.  . is  the  imaginary 
part  (Sine  term  coefficient) . The  j represents  the  first  digit  and  i repre- 
sents the  second  digit  in  each  two  digit  group  of  Figure  3.  The  periodicity 
of  the  Fourier  Transform  allows  duplication  of  the  FFT  output  in  quadrant  I, 


4.  Capt.  Sikula,  G.J. , "Spectral  Signatures  of  Several  Cloud  Types  and 
Information  from  Very  High  Resolution  Visual  Satellite  Radiances  - 
Preliminary  Results" , Sixth  Conference  of  Aerospace  and  AoTonau^irai 
Meteorology,  12-14  November  1974. 
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into  the  II,  III,  IV  quadrants  of  Figure  3. 

There  is  a factor  that  must  be  considered  in  discrete  Fourier  Transform 
processing  because  it  may  bias  the  power  spectra  that  has  been  annularly 
integrated.  This  factor  is  called  aliasing.  Aliasing  is  the  result  of  high- 
frequency  components  inpersona  ting  low-frequency  components.  Aliasing  is 
minimized  by  ocnputing  the  last  half  of  the  FIT  output  by  symmetric  complex 
conjugates.  The  coefficients  for  frequencies  higher  than  the  linear  foldover 
frequency  (or  linear  Nyquist  frequency  equal  to  f J2,  where  fg  is  the  sampled 
domain)  are  computed  by  symmetric  canplex  conjugates  of  frequencies  lower 
than  the  linear  Nyquist.  Since  the  power  spectrum  is  integrated  on  a two- 
dimensional  plane  from  the  output  of  a two-dimensional  EFT,  it  is  acceptable 
to  integrate  the  spectra  annularly  up  to  what  is  termed  the  two-dimensional 
Nyquist  frequency,  given  by: 


/ 


2 , 

(sample  size  of  first  dimension)  + (sample  size  of  second  dimension)^ 

4 4 


Thus  for  a box  in  which  there  are  8x8  picture  elements  (pixels)  the  two- 
dimensional  Nyquist  frequency  is  approximately  6.  Treating  each  two  digit 
number  in  Figure  3 as  a (j,i)  coordinate,  a distance  from  00  (PQ  q)  can  then 
be  computed.  Then  the  terms  (P.  .)  with  the  same  distance  from  P.  _ are 

Dfi  UfU 

sutmed.  Power  terms  that  are  the  same  distance  from  PQ  Q can  be  seen  in  the 
alternating  gray  and  white  bands  of  Figure  4.  Power  terms  that  are  above  the 
two-dimensional  Nyquist  (distance  is  greater  than  6)  are  folded  back  over 
into  lower  bands  by  subtracting  each  digit  from  the  order  of  the  imagery 
array.  That  is  why  the  power  terms  above  Pg  g in  Figure  4 are  in  a band  of 
smaller  distance  from  PQ  Q.  Integration  of  the  FIT  coefficients  in  each 
annular  band  using  this  technique  is  done  by  running  through  the  FIT  output 
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Fig.  3 Schematic  of  Fast  Fourier  Transform  (FFT)  output  power  terms  of 
quadrant  I duplicated  in  quadrants  II,  III,  and  IV.  Each  two 
digit  group  represents  a tvo-dimensional  FFT  coefficient  of 
power  Pj  ^ where  j represents  the  first  digit  and  i represents 

the  second  digit  of  a 25  x 25  nm  box. 
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array  only  once.  Therefore,  duplication  of  the  FPT  output  arrays  as  seen  in 
Figures  3 and  4 was  not  required. 

The  actual  confutation  that  is  made  is  not  power,  but  what  my  be 
termed  normalized  average  amplitudes  (NAA) . Cfcrputations  of  NAA  were 
achieved  by  dividing  the  sunned  power  in  each  of  the  annular  bands  by  the 
number  of  terms  entered  into  the  sum  for  each  band.  The  square  root  of  the 
quotient  is  then  taken  and  normalization  is  achieved  by  dividing  the  average 
amplitudes  by  the  lumber  of  pixels  entered  into  the  FPT.  Normalized  average 
amplitudes  (NAA)  can  be  expressed  by: 

NM(£)  - -vr  /-wrer  pi  <«.*<£> 

Where  P j(f)  ,i(f)  are  Power  spectra  that  are  dependent  upon  the  annular  band 
on  which  they  are  located;  NT(f)  is  the  number  of  P.  ...  . ...  terms  entered 

j \£/  r l (£/ 

into  the  sum  for  each  annular  band;  f is  the  number  of  the  annular  band; 

J and  I are  the  dimensional  orders  of  the  imagery  array.  For  8x8  pixels 
where  J = 8 and  1=8  normalization  is  achieved  by  dividing  by  the  product 
of  J*I  = 64.  Each  of  the  annular  bands  seen  in  Figure  4 has  a unique 
number  attached  to  it  called  the  wave  number  or  frequency.  This  unique 
number  is  actually  the  twD-dimensional  distance  from  the  Pn  term.  In 
essence,  Pq^q  is  the  first  band,  consisting  of  only  one  term,  and  is  said  to 
be  of  frequency  0.  The  first  white  annular  band,  in  Figure  4,  consisting  of 
eight  terms,  is  frequency  1.  The  gray  annular  band  following  the  frequency 
1 white  band  consists  of  twelve  terms  and  is  frequency  2.  This  numbering 
scheme  continues  out  to  frequency  6,  the  outer  gray  band  in  Figure  4,  which 
consists  of  thirty-two  terms. 


6.0  A Study  of  Two  Cloud  Areas 

NAA's  ware  computed  fran  the  coefficients  of  FFT's  of  the  3 nm  resolu- 
tion satellite  imagery  for  20  October  1976.  Two  well  separated  areas,  each 
of  which  contained  a roughly  homogeneous  cloud  field  over  at  least  a 200  rm 
square,  were  selected  for  the  data  sample.  Identification  of  the  clouds  in 
these  areas  was  made  solely  fran  the  only  available  hard  copy  - a small, 
poor  quality,  transparency  of  a DMSP  global  mosaic.  Area  1 was  classified 
as  a mixture  of  stratus  (St)  and  stratocumulus  (Sc)  clouds  and  area  2 was 
classified  as  a uniform  deck  of  St  or  fog.  FFT's  were  performed  on  each 
200  nm  square  area  as  well  as  on  the  64  - 25  nm  sq  (8x8  pixels)  boxes 
that  made  up  the  larger  areas. 

Figures  5 and  6 are  the  means  of  the  NAA  for  64  - 25  x 25  nm  boxes  as 
a function  of  frequency  for  area  1 and  area  2 respectively.  There  are  no 
great  differences  between  the  NAA's  of  the  two  areas,  other  than  at 
frequency  0.  The  difference  in  frequency  0 indicates  that  the  average 
brightness  of  area  1 is  greater  than  area  2.  Figures  7 and  8 are  graphs 
of  NAA  for  each  200  x 200  nm  area  respectively.  Examination  of  Figures  7 
and  8 supports  what  is  seen  in  Figures  5 and  6. 

A statistic  was  computed  to  determine  whether  the  differences  in  mean 
NAA  of  the  two  areas  as  a function  of  frequency  were  significant  using  each 
25  x 25  rm  box.  The  statistic,  T'  measured  u^-  Uj  for  each  frequency  where 
u^  represents  the  mean  NAA's  for  area  1 and  Uj  represents  the  mean  NAA's 
for  area  2.  Table  I lists  the  statistics  T'  as  a function  of  frequency. 

At  the  0.01  level  of  significance  any  T'  larger  than  2.576  or  smaller  than 
-2.576  indicates  that  an  important  difference  in  mean  NAA  exists.  Table  I 
shews  that  there  is  no  significant  difference  in  mean  NAA  for  frequencies  1, 


64  X 64  PIXELS 


FREQUENCY 


FIG.  7 NORMALIZED  AVERAGED  AMPLITUDES  ( NAA)  FOR 
THE  200  x 200  NM  BOX  AS  A FUNCTION  OF 
FREQUENCY  FOR  AREA  I. 


FIG.  8 NORMALIZED  AVERAGED  AMPLITUDES  (NAA)  FOR 
THE  200  x 200  NM  BOX  AS  A FUNCTION  OF 
FREQUENCY  FOR  AREA  2. 
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2,  5 and  6,  however  significant  differences  do  exist  for  frequencies  0,  3 
and  4.  Differences  in  mean  NAA  for  frequencies  3 and  4 iirply  that  in  one  of 
the  areas  there  is  a fluctuation  in  brightness  between  pixels  that  does  not 
appear  in  the  other  area. 

Figures  9 and  10  are  graphs  of  the  mean  variance  of  NAA  as  a function 
of  frequency  for  the  64  - 25  x 25  m boxes  of  area  1 and  area  2 respectively. 
The  variation  of  NAA  as  a function  of  frequency  for  area  2 is  what  would 
be  expected  of  a cloud  deck  with  uniform  tops,  i.e.  essentially  identical 
brightness  between  adjacent  pixels.  Examination  of  Figures  9 arri  10  shows 
that  area  1 has  greater  variability  in  NAA  than  area  2.  Both  figures  show 
the  same  trend;  however  the  variances  are  higher  in  area  1,  indicating  a 

cloud  deck  with  non-uniform  tops  and/or  breaks  above  satellite  sensor  resol- 
ution. 

Table  II  lists  the  statistics,  S^,  S*,  and  F where  sj  are  the  variances 
in  NAA  for  area  1,  are  the  variances  in  NAA  for  area  2,  arri  F is  the 
ratio  of  area  1 variance  to  area  2 variance  as  a function  of  frequency.  At 
the  0.01  level  of  significance  any  F larger  than  1.84  or  snaller  than  0.54 
indicates  that  the  variances  in  NAA  for  area  1 and  area  2 are  very  differ- 
ent for  frequencies  0,  4,  5 and  6.  The  difference  in  variance  for  frequency 
0 means  that  some  of  the  64  - 25  x 25  rm  boxes  in  area  1 have  greater  varia- 
bility in  cloud  brightness  from  one  pixel  to  another  while  area  2 does  not 
exhibit  this  variable  cloud  brightness.  The  ratios  (F)  of  greatest  interest 
and  implication  in  Table  II  are  for  frequencies  4,  5 and  6.  The  reason  is 
that  at  higher  frequencies  the  sampling  rate  of  the  FET  is  at  a higher  cycle 
than  the  lower  frequencies  which  means  that  the  finer  variable  brightness 
in  the  imagery  is  being  sampled.  Thus  if  one  area  has  a finer  variable 
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FIG.  9 VARIANCE  OF  NORMALIZED  AVERAGED 
AMPLITUDES  (NAA)  FOR  64-25x25NM 
BOXES  AS  A FUNCTION  OF  FREQUENCY 
FOR  AREA  I. 
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FIG.  10  VARIANCE  OF  NORMALIZED  AVERAGED 
AMPLITUDES  (NAA)  FOR  64-25  *25  NM 
BOXES  AS  A FUNCTION  OF  FREQUENCY 
FOR  AREA  2. 
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brightness  in  its  imagery  than  another  area  then  the  variance  of  NAA  as  a 
function  of  frequency  of  separate  areas  will  most  certainly  differ.  This 
is  precisely  what  is  seen  in  Table  II  far  frequencies  4,  5,  6,  indicating 
that  the  imagery  for  area  1 is  finer  in  texture  than  area  2.  This  finer 
texture  may  be  the  basis  for  cloud  discrimination  and  thus  the  differences 
in  variances  in  NAA  may  be  a spectral  signature.  In  this  initial  study  it 
would  be  conjecture  to  state  definitively  that  it  is  indeed  a spectral  signa 
ture.  However,  the  important  factor  is  that  a difference  in  imagery  has 
been  detected  in  FFT  processing  of  the  cloud  fields. 


.a,  • 


7.0  Conclusions 


llus  preliminary  investigation  of  using  second  order  statistics  to 
identify  cloud  types  from  1/3  ran  resolution  data  was  based  on  3 m resolution 
data  for  test  purposes.  Nevertheless,  there  is  an  indication  that  processing 
of  Past  Fourier  Transform  coefficients,  as  described  in  Section  5.0,  will 
allow  classification  of  cloud  types  through  a discriminant  function.  As  in 
any  spectral  analysis,  there  is  no  substitute  for  adequate  spectral  and 
spatial  resolution.  For  example,  satellite  data  of  1/3  rm  resolution  would 
undoubtedly  enhance  the  high  frequency  distinctions  of  cloud  distri- 

bution. This  would  give  greater  discriminating  power  among  the  class 
conditional  covariances  of  the  spectra.  Then,  using  a diagonal  covariance 
matrix,  quadratic  discriminants  would  classify  the  imagery  field  into  its 
proper  cloud  category.  This  aspect  of  the  investigation  is  now  under  way. 
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8.0  Future  Work 

This  report  is  a preliminary  study  based  on  the  earliest  data  available 
fran  AFGWC.  Although  3 ran  imagery  was  used,  diffemeces  in  the  cloud  fields 
were  detected  by  processing  FFT  output.  In  the  light  of  this,  it  is 
expected  that  FFT  output  from  1/3  ran  imagery  will  yield  sufficient  differ- 
ences between  cloud  types  to  allow  automated  classification  using  discrimin- 
ant analysis.  Questions  of  how  many  cloud  types  can  be  identified  using 
discriminant  analysis,  their  error  bounds,  and  additional  processing  of  FFT 
output  (such  as  estimation  of  the  class  conditional  covariance  matrices  of 
NAA)  that  will  be  required,  must  remain  open  until  processing  of  1/3  ran 
imagery  is  complete. 

Plans  for  extending  this  work  include  using  coincident  Block  5-D 
EMSP  visual  and  infrared  imagery  for  FFT  processing  and  cloud  type  classi- 
fication through  discriminant  analysis . improvements  on  cloud  classifica- 
tion by  meteorologists  and  looking  at  the  nearest  neighbor  classification 
for  persistence  will  be  considered. 
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TABLE  1 

FREQUENCY  T'  (dQ  = 0) 

* 0 12.3865 

1 -1.3084 

-1.6793 

3 -3.2153 

4 -2.2957 

5 0.3690 

6 -0.0855 
The  statistic  T'  for  corpariscn  of  differences  in  means  far  area  #1  and 

#2(64  - 25  x 25  nm  boxes) . The  null  hypothesis  is  H^:  - Uj  = dQ  and 

the  alternative  hypothesis  is  - Uj  f1  dQ  where  dQ=  0.  The  symbols 

ux  and  are  the  mean  normalized  average  amplitudes  for  area  1 and  area 
2 respectively.  The  critical  region  implying  dQ  f 0 is  t0.005<T' <^0.005 
for  a type  1 error.  Critical  value  for  0.005  with  63  degrees  of  freedom 
is  fc0 . 005  = 2.576. 


TABLE  2 


Frequency 

S 2 

51 

S22 

F=S^/S. 

0 

71.599 

37.583 

1.9051 

1 

0.358 

0.617 

0.5802 

2 

0.262 

0.186 

1.4086 

3 

0.196 

0.123 

1.5953 

4 

0.213 

0.053 

4.0198 

5 

0.157 

0.031 

5.0645 

6 

0.150 

0.069 

2.1739 

The  statistic  F in  ocrparison  of  variance  for  area  #1  arri  area  #2  (64- 

2 2 

25  x 25  nm  boxes) . The  null  hypothesis  is  H^:  = o2  and  the  alter- 

• . . 2 2 

native  hypothesis  is  H^:  = o2  . For  area  1 and  area  2 respectively 

2 2 

and  o2  are  the  variances  in  normalized  amplitudes.  The  critical 
region  implying  o^2  / is  fa/2  < f < l/fa/2  for  a type  I error. 
Critical  values  for  f0  005  and  l/fg/005  are  1.84  and  0.54  respectively. 
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